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Photonuclear reactions at energies covering the giant dipole resonance (GDR) region are analyzed 
with an approach based on nuclear photoabsorption followed by the process of competition between 
light particle evaporation and fission for the excited nucleus. The photoabsorption cross section at 
energies covering the GDR region is contributed by both the Lorentz type GDR cross section and the 
quasideuteron cross section. The evaporation-fission process of the compound nucleus is simulated 
in a Monte-Carlo framework. Photofission reaction cross sections are analyzed in a systematic 
manner in the energy range of ~ 10-20 MeV for the actinides 232 Th, 238 U and 237 Np. Photonuclear 
cross sections for the medium-mass nuclei 63 Cu and 64 Zn, for which there are no fission events, are 
also presented. The study reproduces satisfactorily the available experimental data of photofission 
cross sections at GDR energy region and the increasing trend of nuclear fissility with the fissility 
parameter Z 2 /A for the actinides. 
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I. INTRODUCTION 

In recent years the study of photonuclear reactions has 
attracted considerable interest. Availability of intense 
neutron rich ion beams [l[ will open new vistas in the 
study of nuclei very far away from the valley of stability 
especially in the vicinity of 78 Ni. The use of the energetic 
electrons is a promising tool to get intense neutron-rich 
ion beams. Photofission of Uranium is a very power- 
ful mechanism to produce such radioactive beams (RIB). 
Although the photofission cross section at giant dipole 
resonance (GDR) energy for 238 TJ is about an order of 
magnitude lower than for the 40 MeV neutron induced 
fission, still it is advantageous because the electrons/7- 
photons conversion efficiency is much more significant 
than that for the deuterons/neutrons. In the photofis- 
sion method for the neutron-rich radioactive ion beam 
production, nuclei are excited by photons covering the 
peak of the GDR where the energetic beam of incident 
electrons of ^50 MeV can be slowed down in a tungsten 
(W) converter or directly in the target (U) itself, gener- 
ating Bremsstrahlung 7-rays |2| which can induce fission. 

The aim of the present work is to obtain the pho- 
tonuclear reaction cross sections at energies covering the 
GDR energy region in a framework of a two step pro- 
cess of the nuclear photoabsorption followed by the pro- 
cess of the competition between the light particle evap- 
oration and fission for the excited nucleus. This model 
was used previously for the photonuclear reactions in the 
quasi deuteron (QD) Q energy region successfully 0, [|[ 
whereas in the present work we explore the nuclei excited 
by photons in the GDR energy region which is particu- 
larly important in relation to the production of neutron 
rich nuclei. 



II. THE GDR PHOTOABSORPTION, THE 
NUCLEAR EXCITATION AND FISSION 

In the total photoabsorption cross section at en- 
ergies covering the GDR region, both the Lorentz type 
GDR cross section cr^ DR and the QD cross section <jQ D 
(albeit small) contribute and therefore 
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In the hydrodynamic theory of photonuclear reactions, 
the giant dipole resonance consists of Lorentz line for 
spherical nuclei [H, 0], corresponding to the absorption 
of photons which induce oscillations of the neutron and 
proton fluids in the nucleus against each other, and the 
superposition of two such lines for statically deformed 
spheroidal nuclei [8, corresponding to oscillations 
along each of the nondegenerate axes of the spheroid. 
The lower-energy line corresponds to oscillations along 
the longer axis and the higher-energy line along the 
shorter, since the absorption frequency decreases with 
increasing nuclear dimensions. Therefore, the semiclassi- 
cal theory of the interaction between photons and nuclei 
entails that the shape of a fundamental resonance in the 
absorption cross section is given by 0, @ 
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where Oj , Ei and Ti are the peak cross section, resonance 
energy and full width at half maximum, respectively. 

The QD model [3] is based on the assumption that in- 
cident photon is absorbed by a correlated n-p pair inside 
the nucleus, leaving the remaining nucleons as spectators 
and is proportional to available number of n-p pairs in- 
side nucleus and to the free deuteron photodisintcgration 
cross section o-d(E 7 ). Thus 
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FIG. 1: The plots of cross sections <Jer H as a function of mass 
number Ar of the evaporation residues for 63 Cu and 64 Zn at 
E y =15 MeV, 20 MeV. For medium mass nuclei such as 63 Cu 
or 64 Zn, there are no fission events at these energies. 



fects in the quasideuteron regime of hard photon absorp- 
tion. The free deuteron photodisintegration cross sec- 
tion is given by ad(E 7 ) = 6L2 ^m^ — mb where 

5=2.224 MeV is the binding energy of the deuteron. 
The quasideuteron model Q of nuclear photoabsorption 
is used together with modern rms radius data to obtain 
Levinger's constant L = 6.8-11.2^" 2 / 3 +5.7^- 4/3 of nu- 
clei throughout the Periodic Table and is in good agree- 
ment p"3| with those obtained from experimentally mea- 
sured <jQ d values. 

The recoiling nucleus can be viewed as a compound 
nucleus having the same composition as the target nu- 
cleus but with an excitation energy of E* given by 
E* = m c 2 [(l+2£ 7 /TO c 2 ) 1/2 -l] where E 1 is the 
photon energy and mo is the rest mass of the nucleus be- 
fore photon absorption. This excited compound nucleus 
then undergoes successive evaporation of neutrons and 
light particles or fission. Hence the photonuclear reaction 
cross section a* is a product of the nuclear photoabsorp- 
tion cross section and the statistical decay probability 
and is, therefore, given by = <r;f .^jf where T r and 
r are the partial and the total reaction widths respec- 
tively, and depending upon the case X stands for the 
GDR or the QD contributions. 



TABLE I: Comparison of the calculated photofission cross 
sections for GDR and QD contributions. 
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nuclei 


Fissility 

z 2 

A 


E-y 
MeV 


„GDR 

mb 


„GDR 
f 

mb 


a QD 
w a 

mb 


QD 

a f 
mb 


237 Np 


36.49 


10 
15 
20 


244.0 
392.6 
92.5 


237.0 
391.7 
91.9 


1.165 
5.442 
10.306 


1.131 
5.427 
10.235 


238 yj 


35.56 


10 
15 
20 


269.1 
360.0 
73.7 


186.4 
336.7 
64.9 


1.140 
5.361 
10.189 


0.787 
5.011 
8.967 


232 


34.91 


10 
15 
20 


247.5 
355.6 
75.3 


17.7 
66.5 
15.7 


1.259 
5.680 
10.574 


0.092 
1.080 
2.212 



<r$ D (Ej = J NZa d (E 1 )e- D ^ (3) 

where N, Z and A are the neutron, proton and mass 
numbers respectively, L/A factor represents the fraction 
of correlated n-p pairs and the function e~ D l E ~i accounts 
for the reduction of the n-p phase space due to the Pauli 
exclusion principle. A systematic study of total nuclear 
photoabsorption cross section data in the intermediate 
energy range shows that D = 0.72A 81 MeV [13 which 
for photon energies upto the pion threshold, agrees rea- 
sonably well with the approach based upon phase space 
considerations using Fermi gas state densities that 
conserve linear momentum for the Pauli blocking ef- 



III. COMPETITION BETWEEN LIGHT 
PARTICLE EVAPORATION AND FISSION 

The statistical approach for nucleon and light-particle 
evaporation and nuclear fission is an appropriate scheme 
for calculation of the relative probabilities of different 
decay modes of the compound nucleus. Such statisti- 
cal decay of the compound nucleus is the slow stage of 
the photonuclear reaction. The full Hauser-Feshbach for- 
malism is used in the present calculations, however, for 
pedagogical reasons, we outline the formal calculations 
in the Weisskopf limit only. According to the standard 
Weisskopf evaporation scheme [l5[ , the partial width Tj 
for the evaporation of a particle j = n, p, 2 H, 3 H, 3 He or 
4 He is given by 

^ ^ S + 7^ r ^ <4n V (E) Pj (E*- Bj -E)EdE 
Tr 2 h pcn(E*) J Vj 

(4) 

where Sj, [ij, Vj and Bj are the spin, reduced mass, 
Coulomb barrier and the binding energy of the parti- 
cle j, respectively. cr 3 inv (E) is the cross section for the 
inverse reaction which means the capture reaction cross 
section of the particle j to create the compound nucleus. 
Pcn and pj are the nuclear level densities for the ini- 
tial and final (after the emission of the particle j) nuclei, 
respectively. 

The Bohr-Wheeler statistical approach [l6[ is used to 
calculate the fission width of the excited compound nu- 
cleus. This width is proportional to the nuclear level 
density pf at the fission saddle point: 
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where Bf is the fission barrier height. The diffuse surface 
nucleus Sierk's [l7j fission barriers (Bf) are used for these 
calculations. The decay of the excited compound nucleus 
is treated in full Hauser-Feshbach formalism [l8[ and sim- 
ulated using the Monte-Carlo method [l9j . Finally, in 
order to calculate the fission probability, the total num- 
ber of fission events in a computer run is counted and 
divided by the total number of photoabsorption events. 
Evaporation from excited fission fragments is also taken 
into account. 

The quantitative description of the process is based 
on the liquid drop model (LDM) for nuclear fission by 
Bohr and Wheeler [l6[ and the statistical model of nu- 
clear evaporation originally developed by Weisskopf . 
Accordingly, the probability of fission relative to neu- 
tron emission in the Weisskopf's limit is given by the 
Vandenbosch-Huizenga's equation [2(| 
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where E* = E* - B n and E* f = E* - B f are the nu- 
clear excitation energies after the emission of a neutron 
and after fission, respectively, where B n is the binding 
energy of the emitted neutron. T n and T t are the partial 
widths for the decay of the excited compound nucleus 
via neutron emission and fission, respectively, and the 
parameters a n and a,f are the level density parameters 
for the neutron emission and the fission, respectively and 
Kq = h 2 j2mr\ where m and tq are the neutron mass and 
radius parameter respectively. The emission probability 
of particle k relative to neutron emission according to the 
Weisskopfs statistical model [15[ is given by 
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where is the level density parameter for the emission 
of the particle k, Jk/jn = 1 for k = p, 2 for k = 4 He, 



1 for k = 2 H, 3 for k = 3 H and 2 for k 



3 He. El 



E* — (Bk + Vfc) is the nuclear excitation energy after the 
emission of particle k 21( . B^ are the binding energies of 
the emitted particles and Vk are the Coulomb potentials. 

The evaporation-fission competition of the compound 
nucleus is then described in the framework of a projec- 
tion ang ular momentum coupled [18J Monte- Carlo rou- 
tine [19(. It follows the correct procedure of angular 
momentum coupling at each stage of de-excitation and 
considers statistical emissions of 7, n, p and a particles. 
Any particular reaction channel r is then defined as the 
formation of the compound nucleus via photoabsorption 



and its decay via particle emission or fission. Thus, fis- 
sion is considered as a decay mode. The transmission co- 
efficients for light-particle emissions are determined from 
the optical model potentials. The shift in the coulomb 
barrier during de-excitation is accounted by calculating 
the transmission coefficients at an effective energy deter- 
mined by the shift. The Gilbert-Cameron level density 
22] is adopted in the calculations. The photonuclear 
reaction cross sections a* are calculated using the equa- 
tion = <j*n r /N where n r is the number of events in 
a particular reaction channel r and N is total number of 
events that is the number of the incident photons. 



IV. 



CALCULATION AND RESULTS 



Each calculation is performed with 40000 events using 
a Monte-Carlo technique for the evaporation-fission cal- 
culation. This provides a reasonably good computational 
statistics. The photonuclear reaction cross sections are 
calculated for the GDR and the QD contributions at dif- 
ferent energies for various elements. The Lorentz line 
parameters cr\,E\,Ti and cr 2 , E2, T 2 , for the GDR cross 
sections are taken from Refs.[23|, [24[ and their mean val- 
ues are used for the present calculations. In Table-1, 
the results of these calculations for the GDR and the 
QD contributions for three actinide nuclei are listed and 
compared. In Fig. 1, the cross sections cr%^ R as a func- 
tion of mass number An of the evaporation residues are 
plotted for 63 Cu and 64 Zn at £ 7 = 15 MeV, 20 MeV. For 
medium mass nuclei such as 63 Cu or 64 Zn, there are no 
fission events [25[ at these energies. The variation of 
pGDH an( j a GDR a |. man y gamma-ray energy points over 
the entire GDR bump are detailed in Table-2 and com- 
pared with the corresponding experimental values which 
are extracted from the plots of Ref.[23j. Although the 
calculated <Ja DR values agree quite well with the exper- 
imental values, the photofission cross sections crf DR are 
overestimated implying more neutron evaporation than 
estimated by the calculations. The statistical error and 
standard errors in the Lorentz line parameters [23| can 
not account for the uniform overestimation. On the con- 
trary, at higher energies the situation appears to be re- 
versed and the present calculations slightly overestimated 
the neutron multiplicities Q. However, the present re- 
sults corroborate our earlier findings at higher energies [i[ 
where the fission probabilities for 232 Th, 238 U and 237 Np 
were found to be as high as ~20%, ^87% and ^99% re- 
spectively. The main physics issue probably remains in 
the second stage of the process which suggests perhaps 
the need for more precise experimental measurements in 
the actinide region since little overcounting of the evapo- 
ration neutrons in the measurements can account for this 
discrepancy. 

The upper limits of the reaction cross sections can be 
estimated using aj = a* /N for cases where not a single 
fission event occurred in N events, where N=40000 is the 
number of incident photons. The statistical error Aa^ 
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TABLE II: Variation of the calculated photofission cross sections of actinides with incident photon energy E- 



target 


Cross 
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e.u 


9.U 


1 n n 
10. U 


11. U 


12.0 


13.0 


1 a n 
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lo.U 


lo.U 


1 ( .0 
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nuclei 
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McV 


McV 


McV 


McV 


MoV 


MoV 


McV 


McV 


McV 


McV 


McV 


237 Np 


°a (Calc.) 


mb 


70.50 


128.86 


244.01 


374.63 


390.37 414.86 


443.94 


392.63 


296.09 


213.54 


156.45 




o-j "(Calc.) 


nil) 


69.70 


126.72 


236.82 


354.73 


373.54 


399.18 


441.63 


391.75 


295.39 


212.51 


155.20 




°"a (Expt.) 


mb 




122 


234 


372 


391 


409 


441 


391 


297 








CT^ "(Expt.) 


mb 




66 


134 


219 


197 


225 


284 


259 


178 






238 jj 


o-„ (Calc.) 


mb 


69.39 


133.20 


269.12 


411.46 


393.40 


412.75 


435.88 


360.02 


253.74 


175.90 


126.37 




^""(Calc.) 


mb 


50.76 


87.61 


186.55 


257.68 


241.51 


329.80 


410.53 


336.27 


224.60 


154.63 


104.90 




(7 a ^rjXpr.^ 
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U'J 


1 9^ 
1Z0 


zoo 


4UO 


391 


409 


A 'i 1 


oOy 


ZDZ 


1 70 






af DR (E X pt.) 


mb 


23 


31 


59 


91 


87 


119 


153 


134 


100 


69 


53 




a° M (Calc.) 


mb 


71.28 


131.69 


247.51 


373.36 


389.29 


409.94 


427.28 


355.55 


253.63 


177.37 


128.18 




af> "(Calc.) 


mb 


0.39 


10.28 


17.99 


34.46 


29.48 


54.65 


132.97 


66.58 


48.12 


34.08 


22.41 




a° M (Expt.) 


mb 




144 


266 


375 


394 


412 


431 


350 


250 








t^ M (Expt.) 


mb 




12 


16 


31 


16 


31 


47 


37 


31 







in the theoretical estimates for the photonuclear reaction 
cross sections a* can also be calculated using the equa- 
tion cr^i Aa* — a* [n r ±y/n^\/N which implies that the 
statistical error in the photonuclear reaction cross section 
Aa? = y/a?a?/N. 

V. SUMMARY AND CONCLUSION 

In summary, the photonuclear reaction cross sections 
are calculated for photon induced reactions in the GDR 
energy region within a Monte-Carlo framework for simu- 
lation of the evaporation-fission competition. The calcu- 
lations are performed assuming 40000 incident photons 
for each calculation which provide reasonably good statis- 
tics for computationally stable results. This model was 



used previously for the photonuclear reactions in the QD 
energy region successfully [3, [|| whereas in the present 
work we explore the nuclei excited by photons at the 
GDR energy region which is particularly important in 
relation to the production of neutron rich nuclei. 

The present calculations provide excellent estimates of 
photoabsorption cross section for the actinides while the 
photofission cross sections are somewhat overestimated 
implying more neutron evaporation than expected. The 
light target nuclei such as 63 Cu, 64 Zn, are not very useful 
for RIB production. The contribution to the photofission 
cross sections from the QD process is very small com- 
pared to that from the GDR. The fission channel for the 
medium mass nuclear targets such as 63 Cu and 64 Zn, is 
analyzed and found to be negligible in our calculations. 
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